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NOTE  ON  EVAPORATION  IN  POROUS  MEDIA 
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\R  L Meyer* 

Mathematics  Research  Center 
University  of  Wisconsin 
Madison,  WI  53706 

ABSTRACT  .‘""Vac tors  are  discussed  which  govern  evaporation  of  liquid 
in  the  small  capillaries  of  a porous  medium.  Attention  is  directed  to 
sheet-like  aggregates  from  which  the  vapor  can  escape  with  little 
obstruction.  Marked  temperature  gradients  are  then  found  to  be  confined 
to  close  neighborhoods  of  the  menisci  and  evaporation  is  shown  to  proceed 
in  statiscally  quite  unstable  configurations  under  a dynamic  balance  of 
surface  tension,  local  evaporation  rate  and  viscous  shear.  Estimates 
of  evaporation  rates  and  fluid  velocities  are  given.  The  results 
discourage  constitutive  theories  for  porous  media’ because  mere  size  of 
capillaries,  independently  of  shape  and  chemistry,  isj^ound  to  change  the 
physical  processes  underlying  macroscopic  behavior,  (j  * 

I.  INTRODUCTION.  The  following  study  was  prompted  by  recognition 
that  little  is  known  about  the  physics  of  evaporation  in  fabrics  beyond 
the  guess  that  the  rate  of  heat  supply  may  equal  the  rate  of  latent-heat 
expenditure.  Fabrics  come  in  a great  variety  of  very  different  structures 
and  as  a first  step,  a structure  characteristic  of  "typical"  porous  media 
is  here  envisaged  in  which  the  solid  matrix  is  threaded  by  an  irregular 
network  of  interconnecting,  small  capillaries  along  each  of  which  the 
capillary  bore  varies  greatly  over  relatively  small  distances.  The 
immediate  challenge  is  then  to  isolate  some  of  the  many  interacting, 
physical  processes  in  a single  capillary  in  order  to  distinguish  those 
which  really  govern  evaporation  there;  macroscopic  descriptions  must 
needs  reflect  the  insights  thereby  gained. 

A key  restriction  that  helps  in  dividing  the  difficulties  is  to 
focus  attention  on  sheet-like  media  which  are  thin  in  one  direction,  like 
fabrics,  because  the  escape  of  trie  vapor  is  then  relatively  unobstructed 
and  consideration  of  the  proc. sses  in  the  vapor  can  be  postponed  (to 
Section  VIII).  At  first  sight,  evaporation  might  be  expected  to  be 
controlled  by  the  manner  of  heat  supply,  but  for  capillaries  of  realistically 
small  size,  most  forms  of  heat  supply  have  similar  effects  because  heat 
transfer  across  the  capillary  wall  is  then  always  important.  Since  the 
physical  signposts  diverge,  unless  one  be  quite  specific,  attention  is 
restricted  to  liquids  similar  to  water,  to  pressures  and  temperatures  typical 

• -A  -2 

of  the  outdoors,  and  to  throat  diameters  of  about  10  to  10  cm.  A final 
dividing  step  is  to  start  with  an  unrealistic  configuration  of  geometrical 
and  thermal  symmetry  in  which  only  a single  meniscus  needs  to  be  considered 
(Section  III,  IV). 

Analysis  of  the  simple  thermal  balances  for  that  case  shows  that 
significant  temperature  gradients  can  occur  only  very  close  to  menisci 
(Section  IV),  and  a rough  estimate  of  evaporation  emerges  (Section  IV). 
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It  shows  the  symmetrical  case  to  be  normally  unstable  (Section  V). 

Evaporation  is,  in  fact,  found  to  proceed  in  statically  grossly  unstable 
configurations  under  a dynamic  balance  depending  drastically  on  viscous 
shear.  The  "Haines  Jumps"  in  the  foreground  of  earlier  accounts  [1  - 3] 
can  occur  only  in  much  larger  passages  than  would  appear  realistic  for 
soils,  oil  recovery  or  fabrics.  Instead,  the  normal  state  in  evaporation 
is  one  of  slow  liquid  motion  leaving  the  small  menisci  almost  stationary, 
while  most  of  the  mass-evaporation  occurs  at  the  large  ones  (Section  VI). 

These  results  furnish  a basis  (Section  VII)  for  statistical  estimates  of 
macroscopic  evaporation  rates,  provided  enough  is  known  about  the 
statistical  distribution  of  capillary  throat  sizes.  Such  knowledge 
appears  to  be  an  absolute  prerequisite  for  any  useful  treatment  of  fluid 
motion  in  porous  media  because  viscous  shear,  depends  so  violently  on 
throat  size.  As  a result,  if  two  media  have  the  same  chemistry  and 
identical  shape  for  their  respective  void  passages,  but  differ  in  mere 
geometrical  scale,  then  microscopic  dynamic  balances  can  be  quite  different. 
This  discourages  constitutive  theories  of  porous  media,  of  which  invariance 
to  geometric  scale  is  a basic  premise. 

There  are  other  caveats,  for  instance,  if  the  vapor-air  mixture 
must  pass  through  small  throats,  the  significant  gasdynamical  processes 
must  be  anticipated  to  change  the  evaporation  rate  drastically.  On  the 
other  hand,  there  are  also  many  bits  of  luck,  which  make  a realistic 
fluid  mechancis  of  porous  media  more  accessible.  In  particular,  the 
extreme  magnitudes  of  relevant  combinations  of  physical  parameters  will 
come  to  explain,  by  and  by,  why  errors  by  only  a factor  2,  or  so,  will 
be  treated  so  cavalierly. 

II.  STATIC  EQUILIBRIUM.  When  the  escape  of  the  air-vapor  mixture 
is  unobstructed,  the  pressure  throughout  that  gas  is  effectively  the 
ambient  pressure,  P . Admittedly,  since  surface  tension  promotes 

evaporation,  thermal  equilibrium  between  liquid  and  vapour  requires  a 
gas  pressure  at  their  interface  which  depends  on  the  meniscus  curvature. 
Helmholtz* analysis  [4],  however,  shows  this  to  be  a threshhold  effect,  and 
for  realistic  capillary  bores,  the  thermal  conditions  here  envisaged 
are  well  below  the  threshhold  [5].  With  the  apparent  contact  angle  6 
measured  as  in  Figure  1,  the  pressure  on  the  liquid  side  of  a meniscus 
is  therefore 

T Pa  - 2o/a,  (1) 

where  a denotes  the  local  capillary  radius  and  o sec  ft  the  surface  tension. 

It  will  be  assumed  that  0 < B < n/2,  for  otherwise,  surface  tension  would 
have  kept  the  liquid  out  of  the  porous  matrix.  Body  forces  will  be  neglected. 
Since  it  is  prohibitive  to  take  account  of  all  the  shapes  of  void-passage 
cross-sections  liable  to  occur  in  a porous  matrix,  a will  be  defined  by 

2 

na  = area  of  cross-section;  the  error  in  (1)  is  then  one  of  those  treated 
cavalierly. 

By  (1),  a connected  liquid  column  is  in  static  equilibrium  if,  and  only 
if,  a takes  the  same  value  at  all  the  menisci  bounding  the  column.  The 
equilibrium  is  stable  if,  and  only  if,  a does  not  decrease  with  distance 
along  the  capillary  measured  toward  the  gas  side  at  any  of  those  menisci 
(Fig  2). 


694 


liquid 


K 


Those  are> merely  local  conditions,  however.  All  kinds  of  void- 
passage  shapes  can  be  expected  in  a porous  matrix.  Figure  3 indicates 
one  with  a liquid  column  in  a position  of  stable  equilibrium.  Suppose, 
however,  that  this' liquid  volume  be  reduced  by,  say,  30?o,  e.g.,  by 
evaporation.  The  remaining  liquid  column  cjuld  not  then  find  a stable 
equilibrium  position  anywhere  in  the  passage  segment  shown  : it  must 
move  to  find  a new,  stable  position  elsewhere.  The  need  for  such 
"Haines  jumps"  [1]  in  the  liquid  configurations  has  dominated  the 
literature  on  the  physics  of  fluids  in  porous  media  [1  - 3]  : on 
kinematic  considerations,  they  would  be  expected  to  be  sudden  and 
frequent  so  that  fluid  motion  could  ‘appear  as  a continuous  process  only 
on  a long-term  average  [3],  Analysis  of  the  dynamics  (Section  VI),  however, 
will  show  that  such  notions  can  be  relevant  only  to  unrealistically  large 
passages. 

III.  THERMAL  BALANCE.  Prolonged  evaporation  must  depend  on  an  external 
heat  supply,  and  it  might  be  anticipated  that,  not  only  the  rate,  but  also  the 
manner,  of  that  supply  has  a major  influence  bn  the  evaporation  process. 
However,  since  the  ratio  of  capillary  volume  to  surface  ares  is  proportional 
to  the  capillary  bore,  heat  transfer  across  the  capillary  wall  is  always 
important  in  small  capillaries.  That  transfer  will  be  found  in  section  IV 
to  cause  an  adjustment  of  temperatures  in  the  fluid  and  solid  that  reduces 
the  influence  of  the  manner  of  heat  supply.  To  fix  the  ideas, it  will  be 
envisaged  that  a reservoir  supplies  it  to  the  solid  matrix,  in  the  first 
place,  so  as  to  heat  it  gradually,  from  the  original,  uniform  temperature 
Tg  of  the  whole  medium,  which  may  be  considered  to  be  known,  to  a level 

T1  at  which  it  is  then  maintained.  Specifically,  TQ  may  be  the  outdoor 

temperature,  and  T^  may  be  near  body-temperature,  perhaps  T^  » T^  + 50°F. 

The  properties  of  the  solid  are  outside  the  scope  of  this  study  and  will 
be  assumed  uniform.  This  does  not  imply  a temperature  field  constant  in 
time  and  uniform  in  space,  because  evaporation  makes  the  menisci  act  like 
moving  heat  sinks  of  changing  strength,  but  the  space-average  T of  the 

capillary  wall  temperature  will  change  only  slowly  with  time.  The  stronger, 
local  variation  in  capillary  wall  temperature  can  be  accounted  for 
approximately  by  rough  adjustments  [5]  and  meanwhile,  the  considerations  can 
be  simplified  by  ignoring  the  difference  between  T^  and  the  local  temperature 

of  the  capillary  wall.  That  will  lead  to  estimates  of  evaporation  at 

menisci  in  terms  of  T , whence  estimates  of  the  time  development  of  T by 

w ■ w 

more  global  considerations  will  follow. 

Since  little  can  be  known  about  the  shapes  of  realistic  void-passage 
cross-sections,  the  description  of  evaporation  will  be  simplified  greatly 
by  representing  quantities  in  the  fluid  by  their  averages  over  a capillary 
cross-section  and  ignoring  the  errors  resulting  from  use  of  somewhat 
different  averages  in  different  contexts.  The  thermal  description  is  also 
greatly  simplified,  if  no  liquid  motion  couples  the  processes  at  different 
menisci.  That  is  possible  in  the  idealized  case  of  a capillary  of  radius 
a(x)  even  in  distance  x along  the  capillary,  if  the  temperature  field  is 
similarly  even  and  liquid  fills  a capillary  segment  between  menisci  at 
x = ± £,,  with  a'(£)  > 0 for  static  stability  (Fig  2a).  Such  symmetry  can 
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persist,  at  least  in^principle,  and  will  be  assumed  in  this  Section  and 
the  next  one. 

The  local  thermal  balance  .per  unit  length,  is  then 


2 • 3T* 
7,3  P*C*JT 


_1  ( zx  c j 
' 3x  (lTa  x ax 


) +.-2nah  (T  - T. ) , 

W X/ 


where  the  lefthand  side  represents.the  local  rate  of  increase  of  liquid 
heat  content;,  p and  c denote  density  end  heat  capacity,  respectively, 
and- the  suffix  fi.  iistinguishes  liquid. properties.  The  first  term  on  the 
righthand  side  represents  the  contribution  from  heat  conduction  in  the 
liquid,  and  the  last  term,  that  from  Heat  transfer  across  the  capillary 
wall  on  the  most  rudimentary,  conventional  model  of  a heat  transfer  rate 
per  unit  wall  area  and  unit  temperature  difference  represented  by  a 
constant  transfer  coefficient  h.  With  the  insignificant  fui  .ner 
approximation  of  neglect  of  variations  in  A,  and  c^,  the  balance 

becomes 


3T* 

at 


a_ 

ax 


(a2  21) 
A 3x' 


2h 


P0c 


- (T  - 

Krl 


V 


(2) 


where  h denotes  the  usual  heat  diffusivity,  h = A/  (p^c^).  The  liquid 

therefore  experiences  a typical  heat  conduction  process  with  variable 
effective  diffusivity,  oh  account  of  the  capillary  shape,  and  with  heat 
transfer,  but  without  convection,  on- account  of  the  symmetry. 

A somewhat  different  balance  arises  at  a meniscus.  Since  the  liquid 
and  gas  are  there  envisaged  in  dew-point  equilibrium  to  begin  with,  and 
since  the  gas -pressure  remains  at  the  ambient  ievel  p until  Section  VIII, 

any  heat  reaching  the  meniscus-  will  result  meanwhile  in  evaporation,  but 
not  [2]  iri  a change  of  the  local  temperature  T^  from  its  original  level 

Tg.  Conduction  through  the  liquid,  column  in  x < & contributes  heat  to 

the  meniscus  at  the  rate 

-"a,>3V  > 


if  the  meniscus  is  at  x = A(-t)  and  a denotes  the  capillary  radius  a(£) 

m 

there.  Conduction  through  the  gas  in  x > 11  does  not  contribute  comparably 
because  its  heat  conductivity  is  smaller.  If  the  meniscus  is  markedly 
curved,  heat  reaches  ,it  also  by  direct  transfer  across  a short  segment 
of  the  capillary  wall  and  radial  heat  conduction.  The  wall  area  from  the 
meniscus  contact  line  to  the  position  of  its  apex  (Fig  1)  is  approximately 

2 

2na  with  0 < c??  = secB  - tanS  < 1.  The  rate  of  heat  transfer  across 
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this  area  is 


2ira  a.hM(T  - T ) , 
m 2 M w M 

where  h..  denotes  a value  of  the  transfer  coefficient 
M 

compensate  for  the  error  made  by  confusing  the  local 
at  the  meniscus  with  its  level  further  away. 


h adjusted  [5]  to 
wall  temperature 


Let  S denote  a short  capillary  segment  of  fixed  length  which  is 
stationary  in  a frame  moving  with  the  local,  liquid  velocity  and  which 
contains  the  meniscus  at  present.  The  gas  pressure  and  temperature  are 
constant  in  S,  and  if  Di/Dt  denotes  the  velocity  of  the  meniscus  in  that 
frame,  the  mass-rate  of  evaporation  in  S is 

m = -na  2 p0  D£/Dt  . (3) 

m 

The  mass  of  gas  in  S does  not  increase  at  a significant  rate  because  the 

density  ratio  p /p.  is  about  10”^,  in  the  circumstances  envisaged,  so  that 
g x. 

vapour  leaves  S at  the  same  mass-rate,  and  the  net  rate  of  mass  loss  in 
S is  also  m;  since  no  liquid  enters  S.  By  the  First  Law,  the  net  rate  of 
liquid  enthalpy  loss  in  S equals  ’ue 'rate  of  vapor  enthalpy  loss  from  it 
less  the  rate  of  heat  addition  by  transfer  and  conduction  into  S, 

nam2[2a2hM(Tw  - y - AST^x]  = mL,  (4) 

where  L is  the  latent  heat  per  unit  mass. 
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IV.  THERMAL  LAYER.  The  use  of  these  balances  requires  a nondimensional 
notation,  and  it  is  not  obvious  whether  a single  length  scale  X can  be 
representative  of  the  temperature  field  throughout  a liquid  column.  Even 
without  attention  to  cross-sectional  shape,  the  capillary  is  described 
by  the  two  functions  a(x)  and  a'(x)/a(x)  of  normally  quite  different 
magnitued  and  each  of  which  may  vary  by  orders  of  magnitude  along  a 
capillary.  To  avoid  confusion,  let  attention  be  confined  first  to  a 
liquid  column  segment  adjacent  to  the  initial  meniscus  position  x = i,(0) 
and  short  enough  to  occupy  only  a capillary  segment  characterized  by  a 
single  triplet  of  scales  a^  of  the  capillary  radius,  G,  of  a(x)/a'(x)  and  X, 

of  the  unknown  temperature  variation.  But,  there  is  another  thermal  length  scale, 


n 


8 

3 


£ 

i 


A = (iAaQ/h) 


>2  /v2 


of  decisive  significance  because  A /X  represents  the  ratio  of  the  thermal 
diffusion  scale  to  the  transfer  scale.  How  is  X related  to  the  other 
three  length  scales? 

The  present  model  can  give  no  information  on  how  A compares  with  a^ 

and  G,  but  a more  detailed  calculation  [5]  shows  that  A must  be  anticipated 
to  be  of  the  order  of  the  capillary  radius,  so  that 

A/aQ  = CX/(2aQh) j-  = yh 

is  a parameter  of  order  unity;  a rough  estimate  [5]  is  88  Accordingly, 
X = Sq,  unless  this  be  found  to  imply  that  T^(x)  can  vary  only  on  a longer 
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scale.  The  natural  temperature-difference  scale  is  Tw  - = A,  and  if 

Tw  - T£(x,t)  = T(x,t)  A , 

t is  measured  in  units  of  a time  scale  x,  a/a'  in  units  of  G,  and  x, 
a and  £,  in  units  of  X = a^,  then  the  nondimensional  form  of  (2)  and 


(4)  is 


a0  3T  2a'  ao  3T  32T  1 . , 

nr  3t  a IT  3x  ‘ y^2a 


II.X.-3L  § at  x = t(t), 
3x  . 2 £HT  dt 


respectively,  because 
TU.t)  = 1 

in  this  notation  and  because  the  liquid  is  at  rest;  here  a ^ = o^t^/h  anc^ 
e = c^A/L  « 1 

in  the  circumstances  here  envisaged,  eg,  e = 1/20  for  water  and  A = 50°F. 

Two  different  time  scales  have  emerged  from  the  balances-.  The  shorter, 

2 

aQ  /h,  characterizes  transients  arising  from  imbalance  of  heat  transfer  and 

conduction  which  might  be  anticipated,  eg,  in  a Haines  jump.  It  is  normally 
a rather  small  fraction  of  a second,  eg,  if  h = 10"^cm2/Qec  and  a^  = 10~^cm, 

2 -1 

then  a^  /h  = 10  sec.  The  motion  of  the  meniscus,  on  the  other  hand,  is  on 
the  longer  time  scale 
aQ2/(xe)  . 

The  first  question  must  be  whether  a relatively  stable  evaporation  process 
is  possible,  and  to  examine  this,  the  time  scale  x must  be  identified  with 

2 

3q  /(mc).  If  also  Sg/G  <<  1 , as  one  would  usually  expect,  then  the  lefthand 
side  of  (5)  becomes  unimportant  by  comparison’to  the  righthand  terms,  in 
which  a » 1,  to  the  same  approximation.  The  balances  then  imply 

T = exp  , l(t>  = 1(0)  - lyX  t (6) 

h \ 

[except  when  the  capillary  segment  under  scrutiny  contains  x = 0,  in  which 
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case 


T _ 2 ^(t)/Yh  cosh(x/Y^)]. 

This-  solution  describes  a very  short:  meniscus  layer,  of  thickness  equal  to 
the  thermal  length  scale  A = in  which  virtually  the  whole  process  of 

heat  transfer  to,  and  heat  conduction  in,  the  liquid  takes  place. 

This  conclusion  destroys  the  analysis  sketched  so  far  because  one  of 
its  main  premises  — that  thermal  balances  can  be  formulated  in  terms  of 
cross-sectional  averages  — cannot  apply  to  precisely  the  short  capillary 
segment  containing  the  curved  meniscus  and  all  significant  temperature 
gradients!  Any  tenable  analysis  of  the  temperature  field  must  account 
for  the  geometry  of  the  meniscus,  but  that  depends  mainly  on  matters 
accessible  only  to  vague  speculation,  at  Oest,  namely  the  shape  of  the 
capillary  cross-section  and  the  apparent  contact  angle. 


If  a tenable  analysis  could  be  performed,  n the  other  hand,  it 
would  necessarily  lead  to  the  same  dimensional  groups  and  would  therefore 
also  predict  a dimensional  meniscus  velocity  of  the  form 


A T -T.. 

Y — yii- 


orr 


and  a dimensional  evaporation  time  of  the  form 


y'Vq/^) 

for  a capillary  segment  of  length  a^  What  has  no  rational  support  is  the 
2 

value  (y^+  aj)/Yh  the  nondimensional  coefficient  y predicted  by  (6). 

Thought  about  extreme  cases  indicates  however,  that  the  correct  value  of 
Y cannot  plausibly  be  far  from  order  unity  and  indeed,  that  y = 4 or  5 

cannot  usually  be  very  wrong  To  fix  the  ideas,  therefore,  the  value  y s A 
will  be  adopted  speculatively  for  illustration.  For  water  (with 


m ~ 0.0014  cfn  /sec)  and  z 
give  roughly  the  meniscus 
following  table. 

= 1/20,  various  capillary 
velocites  and  evaporation 

radii  and  lengths  then 
times  listed  in  the 

TABLE  1 

3g(cm) 

'O'2 

io"5 

10"^ 

4eK/aQ  3 x 1Q“2 

3 x 10"1 

3 

cm/sec 

*o  - 

10"1  10"2 

-2  -3 

10  L 10 

cm 

Vo/!teK)  3 

1 1 

7 30 

3 x 10"3  3 x 10"4 

sec 
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In  sum,  the  analysis  has  oeen  wrong  in  everything  but  its  result. 

Most  of  all,  what  has  been  proven,  if  only  by  contradiction,  is  that  the 
temperature  variation  associated. directly  with  evaporation  from  a 
meniscus  must.be  quite  local.  It  follows  that  the  main  results  are  also 
independent  of  some  other  premises  There  can  be  no  significant, direct 
thermal  interaction  between  different  menisci  even  if  the  liquid  moves 

2 

on  the  time  scale  /(he);  the  meniscus  velocity  should  then  be 

interpreted  as  that  of  the  meniscus  relative,  co  the  adjacent  liquid  and 

must  represent  the  scale  of  the  caoillary  radius  a at  the  instanta- 
u m 

neous  meniscus  position.  The  manner  of  heat  supply,  moreover,  can  have 
less  direct  influence  bn  the  local  process  than  might  have  been  thought 
at  first  : apart  from  the  local  dip  in  solid  temperature  at  the  meniscus 
(accounted  for  by  a correct  value  of  y)  the  solid  temperature  background 
of  a capillary  must  reflect  the  macroscopic  scale  of  the  solid  matrix 
as  a whole  and  must  therefore  appear  effectively  uniform  on  the  length 
scale  ag  of  the  local  evaporation  process. 

In  the  first  place,  (6)  applies  only  to  a capillary  segment  in.;  which 
the  radius  differs  by  less  than  an  order  of  magnitude  from  that  at  the 
initial  meniscus  position.  Once  evaporation  has  cleared  that  segment, 
however,  an  analogous  calculation  with  a different  scale  a^  applies  to  the 

next  segment.  Since  narrow  segments  are  seen  to  clear  in  a much  shorter 

time  than  wide  ones,  it  does  not  appear  worth  entering  here  upon  the 

refinement  of  replacing  a from  the  start  by  the  capillary  radius  a (t) 

u m 

at  the  meniscus. 

V.  INSTABILITY.  Before  evaporation,  all  menisci  bounding  a 
connected  liquid  column  must  be  of  the  same  size  (Section  II),  but  in  a 
realistic  porous  medium  a'(x)  cannot  also  be  expected  to  have  the  same 
value  at  different  such  menisci.  The  meniscus  velocity  (Section  IV)  then 
takes  the  same  value  at  all  the  initial  menisci  positions,  but  if  they 
all  started  to  move  with  it,  static  equilibrium  would  be  lost  promptly. 
Surface  tension  acts  towards  restoring  it,  but  the  differences  in  meniscus 
velocity  relative  to  the  liquid  act  in  the  opposite  sense.  Stability  of 
evaporation  therefore  poses  a question  different  from  that  of  static 
stability  (Section  II). 

To  examine  it,  consider  a liquid  column  bounded  by  two  menisci  at 
which  the  gas  pressure,  p , and  temperature,  T = T , remain  the  same, 
but  at  which  the  capillar^  radii  differ.  To  illuminate  the  distinction 
from  static  stability,  suppose  a'(x)  is  monotone  over  the  whole  liquid- 
filled  capillary  segment,  which  contains  a throat  (Fig  4).  Denote  the 
meniscus  positions  by  x = and  x = & < and  the  capillary  radii 

there,  by  a(£+)  = a+  and  a(S._)  = a_  < a+,  respectively;  ie,  the  smaller 

meniscus  is  at  the  lefthand  end  of  the  liquid  column.  By  (1),  surface 
tension  generates  a pressure  difference 


pu+)  - pU_)  = 2a 


driving  the  liquid  towards  the  left. 

If  a'(x)  is  not  too  large,  the  viscous  shear  generated  by  the 
ensuing  liquid  motion  sets  up  a pressure  gradient  related  approximately 
by  Poiseu.ille's  formula 

Q = - •!£  a4  & 

8p  dx 

to  the  mass-flow  rate  Q (counted  towards  the  right),  which  is  independent 
of  x,  by  mass  conservation.  Since  liquid  density  variation  is 
insignificant, 


£ 

ptt  ) - PU  ) = - f + [a(x)]"4  dx, 

+ “ »P  ) 0 


and  since  the  main  contribution  to  this  integral  arises  from  the  throat 
regioh,  it  promotes  clarity  to  write  the  integral  as  in  terms  of 

the  throat  radius  a^.  and  a "throat  length"  fcj..  The  mass-flow  rate 

generated  by  surface  tension  is  then 


_ 3 , 
xpd  t 


a - a 
+ 


4p  a+a_ 


and  the  corresponding  cross-sectional  averages  of  liquid  velocity  are 

2 2 
Q/(npa+  ) at  the  right  meniscus  and  Q/(npa  ),  at  the  left  one. 

Evaporation,  on  the  other  hand,  retracts  the  menisci  into  the  liquid 
with  velocities  YHe/a+  and  yxe/a  , respectively  (Section  IV).  The  center 

of  the  liquid  column  therefore  shifts  at  the  rate 

2 2 2 2 

, Z + Z a—  a r / a,  \ a + a 1 

where 


Ce  = €Hp/(ya^) 

is  a capillary  number  based  on  evaporation  velocity  and  is  normally  very 
small,  if  the  apparent  contact  angle  is  not  close  to  90°  (Fig  1);  for 

-8  -1 

water,  eg,  4agCe  ~ 10  cm.  The  factor  of  Ce  in  the  last  bracket, 
however,  tends  to  be  even  smaller,  as  long  as  a+  and  a remain  large 
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compared  to  the  throat  radius  a^,.  During  such  a phase  of  evaporation, 

the  liquid  column  therefore  shifts  towards  the  right,  ie,  in  the 
direction  opposite  to  that  suggested  by  surface  tension  alone.  Any 
statically  stable  liquid  configuration  with  am»  a^  is  therefore 

unstable  under  evaporation. 

If  the  menisci  were  found  close  to  the  throat,  on  the  other  hand, 
in  the  last  stage  of  evaporation,  then  the  second  term  in  the  last 
bracket  would  be  large,  the  liquid  column  would  move  leftward,  and 
the  smaller  meniscus  would  move  away  from  the  throat. 

In  sum,  most  of  the  evaporation  must  be  anticipated  to  occur  in 
liquid  configurations  that  are  not  static  equilibria,  and  there  might 
te  preferred  positions  for  the  smallest  menisci. 

VI.  DYNAMIC  BALANCE.  For  a realistic  impression  of  evaporation  in 
porous  media,  one  must  therefore  consider  liquid  configurations  far 
from  static  equilibrium,  for  instance,  such  as  indicated  in  Figure  5, 
which  envisages  a situation  that  might  be  seen  in  a snapshot  of  a 
"Haines  Jump"  after  evaporation  has  made  one  meniscus  clear  a throat. 

The  disparity  of  the  meniscus  sizes  then  generates  a marked  pressure 
difference  driving  the  liquid  towards  the  smaller  meniscus,  and  an 
unsteady  liquid  motion  must  be  anticipated.  There  are  two  very  small 
time  scales,  namely  the  liquid  column  length  divided  by  the  sound  speed 

in  it  and  the  time  scale  a^.  /v  of  viscous  diffusion  of  shear  from  the 

-4 

capillary  wall  in  the  throat  region  (Fig  5),  which  is  about  10  sec  in 

water,  if  a.  ~ 10“^cm.  The  evaporation  time  scales  (Section  IV)  are 
much  longer,  and  the  full  viscous  shear  must  therefore  be  expected  to 
have  been  established,  particularly  in  the  throat  region.  The  drastic 
degree  to  which  this  viscous  shear  in  small  capillary  throats  will  be 
seen  presently  to  control  evaporation  illustrates  the  reason  for  the 
prominence  of  Darcy's  law  in  porous  fluid  mechanics. 

The  pressure  imbalance  drives  a mass-flow  rate  Q(t)  < 0,  since  it 
is  directed  towards  the  smaller  meniscus  (Fio  ').  At  the  same  time, 
the  menisci  retract  into  the  liquid  with  the.,  respective  evaporation 
velocities  yHe/a+  and  yxc/a  (Section  IV).  It  will  promote  clarity, 

and  help  to  distinguish  the  more  generic  case  from  that  discussed  in 
the  preceding  Section,  to  exploit  the  disparity  in  meniscus  sizes 
(Fig  5)  to  the  degree  of  neglecting  a against  a . The  larger  meniscus 

2 

is  then  considered  to  move  just  with  the  velocity  dx+/dt  = Q/(npa+  ), 
but  the  smaller,  to  move  with  the  velocity 

dx  /dt  = YH£/a  + Q/(npa  ^)  (7) 

relative  to  a fixed  frame.  The  pressure  difference  is  now  approximated 
as  p+  - p = Zo/'a  and  the  total  shear  stress  is  -8p  QH^/(npa^),  in 
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terms  of  the  "throat  length"  of  Section  V.  In  small  capillaries, 

the  pressure  drop  and  viscous  shear  cannot  come  into  significant 
imbalance,  and  therefore, 

2o/a_  = -8p  Q^/Upa^). 

This  is  the  microscopic  version  of  Darcy’s  law  for  liquid  motion  driven 
through  a small  capillary  by  surface  tension  osecB  at  menisci  of 
disparate  size.  Since  it  relates  Q to  a , (7)  may  be  written 

2 

dx_  r / a 

a-  dt  Y£[  1 4YEHiI^  [ a_ 

and  since  da/dx  < 0 at  the  position  x = x (t)  of  the  smaller  meniscus 

(fig  5),  this  shows  the  approximate  dynamics  to  be  represented  by  an 
equation  of  the  structure 

C1  c!t  = ~^2  ~ 1/°i  = -2Y£a'(x_)  > 0,  (8) 

2 

and  with  increasing  time,  a must  approach 


c2  = oat4/(4YeHpHt)  . 


This  does  not  represent  a strict  equilibrium  because  T , and 

therefore  also  e,  ray  change  slowly  with  time,  and  in  any  case,  the 
liquid  keeps  moving;  but  only  a minor  drift  of  the  smaller  meniscus 
results  therefrom.  In  terms  of  a hybrid  capillary  number 

Ct  = CKp/(oat) 

based  partly  on  evaporation  and  partly,  on  throat  radius,  the  smaller 
meniscus  remains  close  to  the  position  where 

a_/afc  ~ (4YCUt/at)~*  . (9) 


The  capillary  radius  at  the  smaller  meniscus  is  thus  seen  to  depend 

most  of  all  on  the  throat  radius,  indeed,  to  be  proportional  to  , on 

account  of  the  dominance  of  viscous  shear  in  small  capillary  throats. 

For  •-  rough  impression,  for  water,  e = 1/20,  a^  = 10~3cm  and  sec  8 =2, 

the  hybrid  capillary  number  is  Ct  ~ 10"^j  a value  15  of  is  unlikely 
to  be  wrong  by  a large  factor,  and  if  2,^.  = 20at,  then  the  last  formula 
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predicts;  a /a^.  ~ 20.  Figure  5,  accordingly,  gives  a reasonable  impression 

of  s-typical  configuration.  Table  I gives  an.  impression  of  the  liquid 

-3  -2 

velocity  :•  for  a^.  = 10  cm,  eg,  it  is  only  about  10  cm/sec  at  the 
smaller  meniscus,,  and  the  velocity  dx+/dt  of  the  larger  meniscus  relative 
to  a fixed  frame  is  even  less. 

Most  of  the  mass-evaporation,  on  the  other  hand,  occurs  at  the 
larger  meniscus  (Fig  5)  because  its  area  is  larger;  by  (3),  it  is 


m = TtypHea+  = - T^J/L  . (10) 

For  water  and  a+  = 1mm,  eg,  it  would  be  about  10“^gr/sec.  In  turn,  an 

impression  of  the  dependence  of  the  wall-temperature  level  T upon  the 

w 

external  heat  supply  begins  to  emerge,  because  the  rate  of  that 

supply  per  meniscus  is  mL.  In  a liquid  column  bounded  by  just  two 
menisci  of  disparate  size,  the  smaller  expends  relatively  little  of 
this,  so  that  (10)  shows  the  external  heat  supply  to  such  a column 
to  depend  only  on  Tw  and  on  the  capillary  bore  at  the  larger  meniscus. 

The  critical  importance  of  capillary  size  merits  re-emphasis 
here.  The  pressure  difference  due  to  surface  tension  is  proportional 

_1 

to  affl  , and  so  is  the  meniscus  velocity  relative  to  the  liquid,  but 

the  pressure  drop  due  to  viscous  shear  is  proportional  to  a^.  . 

Accordingly,  if  two  porous  samples  be  compared  which  are  identical  in 
regard  to  chemistry  and  to  shape  of  the  void-passages,  but  differ  by  a 
factor  10  in  the  size  of  those  passages,  then  the  dynamic  balances  for 
them  differ  by,  essentially,  a factor  10^  and  therefore,  the  fluid 
physics  in  the  two  samples  may  be  quite  different.  That  contrasts 
strongly  with  constitutive  theories  of  porous-media  mechanics,  of  which 
invariance  under  mere  change  of  geometric  scale  is  a main  premise. 

It  appears  doubtful,  therefore,  that  a substantive  description  of 
fluid  mechanics  in  porous  media  is  obtainable  without  some  insight  into 
the  fluid  dynamics  on  the  microscopic  level,  whence  macroscopic 
behavior  must  spring. 

VII,  MACROSCOPIC  IMPLICATIONS.  For  an  impression  of  global 
evaporation  in  a porous  medium,  a thermodynamically  steady  phase  may 
need  to  be  distinguished  from  an  initial,  transient  phase.  The  later 
phase  is  characterized  by  essential  equality,  at  any  time,  of  the 
rates  of  external  heat  supply  and  of  latent-heat  expenditure.  The 
global  mass-rate  of  evaporation  is  then  immediately  known  and  a lower 
bound  ip  of  total  evaporation  time  can  be  deduced  as  that  which 

evaporation  of  the  initial  liquid  mass  would  need  under  such  conditions. 


Whether  such  a late  phase  occurs  at  all,  or  at  the  other  extreme, 
whether  the  transient  phase  is  of  no  importance,  must  be  judged  from 
comparison  of  tg  with  the  time  scale  of  the  transient  phase.  The  latter 

may  be  one  of  three  scales,  of  which  the  first,  , characterizes  the 

rate  at  which  the  external  reservoir  can  communicate  heat  to  the  porous 
aggregate  and  a second,  x2,  characterizes  the  rate  at  which  heat  can 

be  distributed  through  the  solid  matrix.  Neither  of  these  is  within 
the  scope  of  this  account,  but  the  time  scale  x^  of  liquid  response 

to  the  heat  supplied  to  it  can  be  predicted  on  the  basis  of  the  present 
results,  if  adequate  knowledge  of  the  statistical  distribution  of  void- 
passage  sizes  is  at  hand. 

Indeed,  if  even  a relatively  small  number  of  large  passages  thread 
the  porous  aggregate,  then  all  the  "action"  will  occur  in  them,  and 
Haines  Jumps  may  be  there  observable,  while  the  rest  of  the  medium  remains 
essentially  inert  until  the  by-passes  have  cleared  to  an  extent  making 
them  effectively  a part  of  the  outer  boundary  of  the  medium. 

On  the  other  hand,  if  enough  passage  throats  of  sufficiently  small 
size  are  distributed  sufficiently  well  through  the  aggregate,  then  they 
will  anchor  the  liquid  and  permit  only  creeping  motion.  The  time  scale 
Tj  would  then  be  expected  to  be  essentially  that  of  transition  from 

static  stability  to  dynamic  balance,  which  (8)  shows  to  be 

t3  = c^g/Cne)  = (a0/e)2/[2yx|a'(x_)|]. 

Since  it  is  seen  to  depend  most  of  all  on  ag  and  e,  a useful,  macroscopic 

estimate  of  this  scale  requires  both  a judicious  choice  of  e between 
0 and  its  level  in  the  steady  stage,  and  also  a statistically  valid 


corresponding  measure  3g  of  capillary  radii  at  the  small  menisci  can 

be  deduced  by  (9).  For  a very  rough  impression  (which  may  well  be  mis- 

_2 

leading  in  regard  to  specific  cases),  if  e = 1/50,  ag  = 10  cm  and 

| a ' ( x_ ) | = 10  were  appropriate,  then  for  water,  x^  ~ h min. 

If  t3  should  turn  out  to  be  much  longer  than  ig,  and  x2,  the 

microscopic  dynamics  of  the  main  phase  of  evaporation  would  be  that 
described  in  Section  V.  The  stability  analysis  there  given  is  not  linearized, 
and  its  results  therefore  apply  to  the  whole  transition  from  static 
stability  to  dynamic  balance.  Of  course,  once  the  initial  configuration 
of  static  stability  has  been  left  well  behind,  the  much  simpler  approxima- 
tion of  Section  VI  becomes  adequate.  The  translation  here  of  the  micro- 
scopic description  into  a predictive  algorithm  on  the  macroscopic  scale 
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may  be  -premature , however,  because  its  usefulness  is  likely  to  depend 
critically  on  a more  precise  knowledge  of  the  statistical  distribution 
of  void-passage  sizes  than  appears  available  tg-date  for  any  real  sample. 

vm.  VAPOR  TRANSPORT.  For  evaporation  in  dynamic  balance,  the 
time-dependence  of  the  processes  in  the  air-vapor  mixture  may  also  be 
expected  to  amount  to  no  more  than  a slow  drift  leaving  the  processes 
quasi-steady.  Transfer  from  the  capillary  wall  will  heat  this  gas, 
but  if  the  pressure  differences  in  it  are  insignificant  because  no  small 
throats  obstruct  its  passage,  then  the  attendant  density  change  will 
not  be  worth  accounting  for,  at  the  temperature  levels  here  envisaged. 

2 

Accordingly,  the  volume  flow  rate  ira  u of  gas  will  also  be  considered 
constant  along  the  passage. 

The  mass  evaporated  consists  of  vapor,  but  the  gas  flow  moves  the 
air-vapor  mixture  and  must  therefore  be  accompanied  by  diffusion  of 
vapor  and  air  into  each  other.  The  gas  at  the  meniscus  must  be  at  its 
dew  point,  so  that  the  partial  vapor  pressure  there  is  the  saturation 
pressure  at  the  meniscus  temperature  T^.  For  water  vapor,  e.g.,  that 

partial  pressure  is  about  1/40  (or  3/40)  at  T^  « 293  (or  313)°  K,  and 

the  gas  even  at  the  meniscus  then  consists  almost  entirely  of  air. 

Since  the. partial  vapor  density  .is  even  smaller  [4],  the  diffusion  of 
the  vapor  is  adequately  approximated  by  the  standard,  linear  model  of 
Fick's  law,  Jv  = - 9Vpv,  for  the  vapor  flux.  With  unsteadiness  already 

neglected,  the  same  mass-flow  rate  of  vapor  must  cross  every  capillary 

2 2 

so  that  na  upy  - ua  9dpv/dx  is  independent  of  x.  It 


(C-5J/q 

e e ’ 


where  subscripts  m and  e distinguish  respective  values  at  meniscus  and 
capillary  exit, 

5 = f Ca(s)]"2  ds, 

Jx 

m 

and  q = 9/um  is  a diffusion-length  scale  based  on  the  gas  velocity  um 

at  the  meniscus.  Most  of  the  diffusion  therefore  occurs  within  a £- 
distance  q of  the  exit,  by  contrast  to  the  heat  transfer,  most  of  which 
occurs  fairly  close  to  the  meniscus. 

The  process  can  be  radically  different,  however,  if  the  gas  must 
pass  through  a small  throat.  Let  u denote  again  the  cross-sectional 
average  of  the  gas  velocity  and  let  subscripts  £,  g,  v,  a,  m and  t distin- 
guish reference  to  the  liquid,  gas,  vapor,  air,  meniscus  and  throat, 
respectively.  Then  from  the  estimate  of  meniscus  velocity  relative 


cross-section, 
follows  that 
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to  the  liquid  in  Section  IV,  amum  = yeg^p^/p^,  approximately,  and  since 

this  is  independent  of  meniscus  size,  so  is  the  Reynolds  number  Rem  = 

2 

= a u /v  of  the  gas-flow  at  the  meniscus.  Since  v - 0.'i5  cm  /sec 
m m a a 

and  p„/p  - 10"5  for  water,  Table  1 (Section  IV)  indicates  Rem  = 2 to 
£.  v m 

2 

be  a rather  typical  value.  The  mass- flow  rate  m = Tia  p^u  is  independent 

of  x in  near-steady  evaporation,  and  apart  from  the  influence  of  density 
changes,  the  local  Reynolds  number  Re  = au/v  of  the  gas-flow  varies 

in  proportion  to  a /a(x).  For  most  plausible  values  of  a /a, , the  gas- 
lit m t 


flow  therefore  remains  laminar  even  in  a throat,  and  the  pressure  drop 
can  again  be  estimated  from  Poiseuille's  formula, 

3 2 3 

a dp/dx  = - 8uam/(na  Pg)  = - 8^aamum/a  > 

where  8pgamum  is  independent  of  meniscus  size  and  typically,  - 5 x 10-^  gr 

-7  2 

when  u =2x10“  gr  sec/cm  . 
a 

-2 

If  now  a = 10"  cm,  to  fix  the  ideas,  then  a /a.  = 10  yields  a 
m ’ ’ m t 

value  of  5 x 10-"5  atm  for  |a  dp/dx | at  the  throat,  and  the  pressure 
drop  is  insignificant.  If  a^a^  - 100,  however,  then  the  estimate  suggests 

a value  of  5 atm  for  |a  dp/dx ) at  the  throat,  and  not  only  the  estimate, 
but  clearly  also,  most  of  the  premises  and  assumptions  of  this  Note, 
collapse.  If  the  evaporation  estimates  remained  valid  when  the  gas 
must  pass  through  very  small  throats,  they  would  imply  major  gasdynamical 
effects  in  such  throats,  the  work  expended  on  them  would  play  a major 
role  in  the  thermodynamic  balances,  and  the  meniscus  tem"'>rature  T^ 

could  not  be  expected  to  be  close  to  the  initial,  ambient  temperature 
Tq ; the  physics  of  evaporation  would  be  quite  different  from  that  here 


described.  Accordingly,  the  microscopic  physics  of  evaporation  may 
depend  rather  drastically  on  whether  a porous  medium  is  formed  into 
a sheet-like  or  ball-like  aggregate... 


The  great  sensitivity  of  quantitative  estimates  to  void-passage 
size  suggests  that  a profitable  discussion  of  fluid  mechanics  in  porous 
media  may  need  to  relate  to  quite  specific  circumstances.  In  particular, 
if  thought  returned  to  fabrics,  it  appears  unlikely  that  the  same  micro- 
scopic physics  could  describe  evaporation  from  wool,  gore-tex  or  pile. 
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